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Abstract 

The behavior in an external magnetic field is studied exactly for a wide class of 
multi-chain quantum spin models. It is shown that the magnetic field together with 
the inter-chain couplings cause the commensurate-incommensurate phase transitions 
between the gapless phases in the groundstate. The conformal limit of these models is 
studied and it is shown that the low-lying excitations for the incommensurate phases 
are not independent, because they are governed by the same magnetic field (chemical 
potential for excitations). A scenario for the transition from one to two space dimen- 
sions for the exactly integrable multi-chain quantum spin models is proposed and it is 
shown that the incommensured phases in an external magnetic field disappear in the 
limit of an infinite number of the coupled spin chains. The similarities in the external 
field behavior for the quantum multi-chain spin models and a wide class of quantum 
field theories are discussed. The scaling exponents for the appearence of the gap in the 
spectrum of the low-lying excitations of the quantum multi-chain models due to the 
relevant perturbations of the integrable theories are calculated. 

PACS numbers: 75.10. Jm,11.10.Kk 
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1 Introduction 



There has recently been considerable interest on low-dimensional quantum correlated spin 
and electron systems. These systems, especially one-dimensional (ID), manifest the specific 
features of, e.g., magnetic behavior at low temperatures, which are absent for the standard, 
conventional 3D magnetic systems. Spin systems usually manifest ID behavior for the 
temperatures higher than the temperature of the 3D magnetic ordering, but lower than 
the maximum characteristic energy of the interaction between spins, i.e. in our case the 
intra-chain spin-spin coupling. The origin of such specific features is the enhancement of the 
quantum fluctuations of the ID systems due to the peculiarities of the ID density of states 
together with the quantum nature of spins. 

Moreover, during the last decade a large number of new quasi ID spin compounds were 
created and studied experimentally. These compounds manifest at low temperatures the 
properties of a single or several quantum spin chains weakly coupled to each other EL It 
is strongly believed that this class of compounds will provide the new information on the 
transition from ID to 2D in quantum many-body physics. It is very important, because the 
2D quantum many-body physics has been a challenge for both theorists and experimentalists 
since the beginning of the study of low dimensional quantum systems. On the other hand, 
the advantage of the ID theoretical studies is the possibility of obtaining exact solutions 
by using non-perturbative methods, which are difficult to apply for the higher- dimensional 
quantum many-body models. The results of the exact calculations of the ID models can 
serve as testing grounds for the use of perturbative and numerical methods in more realistic 
situations. 

Recently several exactly solvable models |3], |, [5| have been introduced, in which the 
zigzag-like interaction between two quantum spin chains was studied exactly using the Bethe 
ansatz technique ||. This method is widely known by now, see e.g., the recent monography 
[[/J and references therein. The Bethe ansatz method permits exact calculation of the static 
characteristics of quantum many-body systems, such as the groundstate behavior, the influ- 
ence of an external magnetic field, and the thermodynamic features of e.g., the temperature 
dependencies of the specific heat, magnetic susceptibility, etc. These results should apply 
to more realistic systems, but it is not obvious how the interactions between the chains 
modify the answers. The mean-field like approximations for the inter-chain couplings are 
not sufficient, because the mean field approach in any version already implies the existence 
of the (sometimes hidden) order parameter. It is, unfortunately, also unclear whether the 
numerical calculations, which can be directly applied for the quantum many-body systems of 
very small sizes by now (say, at most several tens of sites) describe well the properties of the 
real systems, in which, even in quasi-lD ones, the number of sites is at least of order of 10 8 
or higher. On the other hand, it must be admitted that some features of the exactly solvable 
ID models are far from what is observed experimentally, but these unrealistic features of the 
ID models are known and simple to recognize. 

The behavior of the multi-chain spin systems in an external magnetic field is especially 
interesting, see e.g., |8], |], [10], |5| because of (i) the possibility of experimental observations 
due to recent progress in the high magnetic field measurements, and (ii) very interesting 
theoreticallly predictable effects which are possible to recognize in experiments, such as phase 
transitions in the external magnetic field. However several important issues are far from being 
solved in the quantum two-chain spin models. For example, there are three questions that 
need to be answered: (1) Are the properties of those exactly solvable two-chain spin models 
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unique or it is possible to say something about the more general class of two-chain quantum 
spin models? (2) How are the multi-chain quantum models connected to the 2D many-body 
systems, i.e. what is the scenario of the transition from ID to 2D when one increases the 
number of coupled chains while keeping the conditions of integrability? and (3) What will 
happen with the behavior of the non-integrable multi-chain spin models if one goes beyond 
the framwork of integrability i.e. adding some perturbations to the exactly solvable model? 
(For example, Ref. implies that namely the spin chirality, which separatly breaks the 
time-reversal and parity symmetries in the two-chain integrable model flUf , is the reason 
for the emergence of the additional phase transitions in an external magnetic field for the 
two-chain spin ~ model as compared to the single-chain system). 

The goal of this paper is to answer these questions. First, we re- visit the exactly integrable 
two-chain spin ~ model and show that the inclucion of the magnetic anisotropy of the "easy- 
plane" type, with which the system stays in the quantum critical region, will not drastically 
change the behavior in an external magnetic field but will shift the critical values of the 
magnetic fields and intra-chain couplings at which the phase transitions occur and will 
affect the critical exponents. We will show that these two-chain spin models share the most 
important features of the behavior in an external field with the wide class of the (1+1) 
quantum field theories. Next, we will introduce the higher-spin realizations of the two-chain 
spin models, e.g., investigating the important class of ID two-chain quantum ferrimagnets 
with different spin values at the sites of each chain. We will also investigate the behavior 
of the multi-chain exactly solvable spin models in an external magnetic field and show how 
the additional phase transitions arising due to the increasing number of chains vanish in the 
quasi-2D limit. Finally, we will show how the relevant deviations from the integrability, e.g., 
the absence of the terms in the Hamiltonian which separately break the parity and time- 
reversal symmetries give rise to gaps in the spectra of low-lying excitations of the multi-chain 
quantum spin systems and we will calculate the scaling exponents for the gaps. 

The paper is organized as follows. In Section 2 we re-visit the exactly solvable two- 
chain uniaxial spin model 0] to remind the reader of the main steps of the Bethe ansatz. 
The investigations of Refs. |], of isotropic spin 1 two-chain models are generalized in 
this section for the case of uniaxial magnetic anisotropy. The calculations in this section 
are rather simple, but we write them in detail because they provide the basis for the more 
nontrivial generalizations of this class of models, and will be used in the following Sections. 
In Section 3 we point out the similarities between the behavior of the uniaxial two-chain 
quantum spin models and a class of quantum field theories (QFT) in an external magnetic 
field, predicting new phases for the QFT. In Section 4 we introduce the SU (2) generalization 
of the integrable two-chain model for higher values of the site spins (possibly different) in each 
chain, i.e., quantum ferrimagnet. We point out the similarities of the quantum ferrimagnet 
with the QFT with a nonzero Wess-Zumino term and predict new phases for the latter in 
an external magnetic field. We derive the integral eguations for the critical exponents. In 
Section 5 we consider the multi-chain quantum spin model and discuss how the external 
field behavior of the integrable multi-chain models is changed when the number of chains 
is increased while preserving the exact solvability. In Section 6 we briefly sketch how the 
deviations from integrability change the magnetic and low-temperature properties of this 
class of multi-chain quantum spin systems. The paper is closed with a discussion of the 
main results and some conclusions. 
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2 Two-chain uniaxial quantum spin model 



A common property of some of the Bethe ansatz solutions is the presence of shifts 9j of the 
spectral parameter A for the associated transfer matrix of an algebraic version of the Bethe 
ansatz (the Quantum Inverse Scattering Method, [QISM] 0). Those shifts also appear in the 
Bethe ansatz equations (BAE) for the quantum numbers called rapidities, which parametrize 
the eigenf unctions and eigenvalues of the Hamiltonians. Hence, the distributions of the 
rapidities are also affected by the shifts. An interesting property is connected with those 
shifts: depending on their values and the external magnetic field, even for (quasi)particles 
of the same type, additional minima may appear in distributions of the rapidities. These 
additional minima also result in the nonmonotonic behavior of the dispersion laws of the 
low-lying excitations. Also, they provide additional Dirac seas for low-lying excitations, 
changing the structures of the physical groundstates of the models. These additional minima 
determine the special behavior of the models in the external magnetic field [[3], |S], ID, In 



particular, the appearence of the new phases and new phase transitions is namely due to the 
emergence of these new minima in the distributions of the quantum numbers. 

To set the stage, let us first remind the reader about the main steps of the QISM. The 
common feature of the Bethe ansatz solvable models is the factorization of a monodromy 
matrix (the ordered product of all two-particle scattering matrices, which depend on some 
spectral parameter) 0. Exact (Bethe ansatz) integrability requires exclusively elastic scat- 
tering between (quasi)particles. For such a theories two-particle scattering matrices and 
//-operators satisfy the Yang-Baxter relation |^2, |7]. In turn, the factorization of the mon- 
odromy matrices garantees that they satisfy the Yang-Baxter equations, too. The transfer 
matrices of the associated statistical problem are traces over some additional, auxiliary sub- 
space, of monodromy matrices f7|. The most important feature of transfer matrices with 
different spectral parameters is their commutativity. The necessary and sufficient condition 
for this is the validity of the Yang-Baxter equations for two-particle scattering matrices and 
hence for monodromy matrices. The commutativity of transfer matrices implies that one can 
construct an infinite number of integrals of motion, which commute with one another and 
with the transfer matrix. Therefore the exact integrability is proved. Usually the structure 
of these integrals of motion is determined by their locality. For instance, the best known of 
series of integrals of motion is the series of derivatives with respect to the spectral parameter 
of the logarithm of a transfer matrix taken at some special value of former [[/]]. Locality 
means that for the first derivative of the logarithm of the transfer matrix (usually called the 
Hamiltonian of the lattice system) only short-range particle-paricle interactions contribute. 

In this paper we will see that namely the aforementioned shifts of the spectral parameters 
yield new phases in the groundstate behavior in an external magnetic field of a wide class of 
exactly solvable models, quantum spin multi-chain models and QFT. We will show that in 
the conformal limit these phases of the lattice models correspond to one Wess-Zumino-Witten 
(WZW) model or to several of them with dressed charges (proportional to the compactifica- 
tion radii) of scalar or matrix types for each of the phases, respectively. 

Let us start with the form of the Bethe ansatz equations (BAE) for the set of rapidities 
{ u a\a=i- m this paper we will concentrate only on the critical, "easy-plane" type of the 
magnetic anisotropy for the antiferromagnetic spin multi-chain models, 0<7<7r/2(7 = 
ir/q, q integer, parametrizing the magnetic anisotropy), and the repulsive interactions in 
QFT. This correspond to hyperbolic or rational solutions of the Yang-Baxter equations 
for two-particle scattering matrices, or to U(l) and SU(2) symmetries of the scattering 
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processes, respectively. For the simplest case of one shift 6, which is connected to the two- 
chain quantum spin models and most of QFT, the BAE have the form (here we use more 
general hyperbolic parametrization first; for the rational limit see below) Q: 

M 

Y[e^(u Q ±6) = e™ M J] e 2 (u a -Up), (1) 

± p=l,pjta 

where N± are the numbers of sites in each of spin chains, e n (x) = sinh(x-M7f )/ sinh(x — Z7§) 
and M is the number of down spins. The shift 9 determines the inter-chain coupling constant 
for two-chain quantum spin ~ models [f|, [11], |13| . Please note that the Bethe ansatz equations 
are just the quantization conditions for the rapidities, which parametrize the eigenwaves and 
eigenvalues of the many-body quantum model. The Hamiltonian is the first derivative of 
the logarithm of the transfer matrix (note that the transfer matrix of the two coupled spin 
chains in this integrable model is the product of two "standard" transfer matrices of each 
chain with the spectral parameters A ± 9 |TTJ ) : 



H\/2 = . - 2 n r^r~ ( sinh 2 9E(S n ^S n +i,i + 5 n 2-^1+1,2) + 

smh 9 + sin 7 „ V 

2 sin 2 7 IS nil (S n>2 + Sn+1,2) + 2sin7sinh6'(J5 ? „ + i i 2 - J S n>1 )[S n+ljl x S n>2 ]j , (2) 

where diag(a, b, c) is 3 x 3 diagonal matrix, 

E = diag(l, 1, cos 7) , 

I = diag(cosh 8, cosh 9, cos 7) , 

J = diag^cosj, cos 7, cosh 9) , 

and [. x .] denotes the vector product. Please note that the sum runs over n to iV + for the 
chain with spins S Hi i and to AL for the chain with spins SV^- The parameter 9 determines 
the intra-chain coupling in our two-chain spin model. For 9 = the Hamiltonian and 
BAE coincide with the ones for the single "easy-plane" antiferromagnetic spin | chain of 
length + AL with the only nearest neighbour interactions in it. The eigenvalue of the 
Hamiltonian (energy) is parametrized as the function of the rapidities as follows: 

M 

E = sm 7 Y,Y, N±{ei{u a ±9)+ e^(u a ±9)) + E , (3) 

± a=l 

where Eq is the energy of the vacuum (ferromagnetic) state (with M = 0). The isotropic 
5't/(2)-symmetric antiferromagnetic quantum spin two-chain model JTXJ, [13], |9|, [1(| can be 
obtained from the uniaxial (t/(l)-symmetric) one of Eqs. (|T])-([J) by the simple change of 
variables in the limit: u a —>■ 'jUa, A — > 7A, 9 — > 7$, 7^0. (The last limit corresponds to the 
rational, 5 , [/(2)-symmetric solution of the Yang-Baxter equations for two-particle scattering 
matrices). The two-chain isotropic (S'f/(2)-symmetric) spin | Hamiltonian obtained in this 
limit from Eq. (|) takes the form H OS £l g, E3 : 



His ~ ~i — 7~no 5^ (@ 2 {S n) lSn+l,l + 5^,25^+1,2) + 2S n t l(S n 2 + S n+ i 2) + 

1 + 9 2 n \ 

29(S n+lt2 - S nA )[S n+1A x S n ,2]) . (4) 
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The summations over n runs to N± for each kind of spins, respectively. Note that for 9 — > oo 
Eqs. (|J) and BAE recover the Hamiltonian and BAE of two decoupled spin | chains of lengths 
N± with the only nearest neighbour interactions in each of the chains. 

The solution to the BAE Eqs. ([I]) is usually obtained in the thermodynamic limit 
(N±,M — > oo, with the ratio M/(N + + AL) fixed). Here instead of the discret set of 
rapidities one introduces the distribution of a continuous density of rapidities. The ground- 
state corresponds to the solutions of the BAE with negative energies, i.e., it is connected 
with the filling up the Dirac sea(s) for the model. For the "easy-plane" antiferromagnetic 
two-chain spin \ model the groundstate corresponds to the filling of the Dirac sea for the real 
rapidities, i.e., no spin boundstates have negative energies. In the thermodynamic limit the 
real roots of Eqs. (Jf]) are distributed continuously over some intervals, which determine the 
Dirac seas of the model. The set of integral equations for the dressed densities of rapidities 
u a (p(u)) and dressed energies of low-lying quasiparticles (s(u)) are (see, e.g., Ref. M for 



the standard procedure of deriving these integral equations from the BAE and Refs. [[LI], [13 
for the isotropic two-chain spin | model): 



p(u) + dvK(u - v)p(v) = £ -fpl (5) 

and 



A^ 4 

dvK{u — v)e{v) = h — } j 
KQ) ± 

where the kernels of integral equations are 



r /v 

e(u) + / dvK(u - v)e(v) = h - £ -^4 , (6) 



= dlne 2 {u) = sin(2 7 ) 
K ) 8u " 27r[cosh( M ) -cos(2 7 )] ' K ) 

and h is an external magnetic field. The values 

d\ne 1 (u±6) _ dp\(u) = sin7 

P±[U) du ~ du 27r[cosh(«±^) -cos( 7 )] [ ) 

are bare densities of the rapidities, and 

o , . sin 2 7 

£±[ u) = h - ______ (9 ) 

are bare energies (here "bare" corresponds to noninteracting particles, and the interaction 
"dresses" them as usual 0). The integrations are performed over the domain (Q), deter- 
mined in such a way that the dressed energies inside these intervals are negative. The limits 
of integrations are determined by the zeros of the dressed energies, and are the Fermi points 
for each sea. The analysis of the integral equations Eqs. in an external magnetic field 

shows that in general, for some values of 9 and h, there can be one Dirac sea (it corresponds 
to one minimum of the bare density of rapidities and, hence to one minimum of the bare 
energy). On the other hand, for higher values of 9 and for some domain of h two Dirac 
seas of the same type of (gapless, see below) excitations are possible (for two minima of 
the bare energies of the rapidities and thus two minima of the bare density). Note that for 
9 — » oo at fixed N± all the roots of the integral BAE separate into two sets of "right-" and 
"left-moving" seas, centered at ±9, respectively. 
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Here we briefly re- visit the analysis of Refs. [§, [H| , but for the case of the uniaxial two- 
chain model. Analytic solutions to Eqs. (§)-(§) can be easily obtained in closed form in the 
limit of zero field and equal lengths of the chains N + = N_. The simplest non-trivial exited 
quasiparticle (spinon) is a hole in the Dirac sea for real rapidities with the quasimomentum 



p(u ) = 2 arctan ( ^A™^ 7 ? ) 

V cosh(7rt77) / 



(10) 



where uq is the spinon's rapidity. Note that due to topological reasons such particles have 

I5[. The energy of this spinon is 



to exist in pairs for the S'[/(2)-symmetric case, etc. (Tj 
given by 

/ v . dpjup) 

e(u ) = - sin 7— . (11) 

It can be rewritten as function of the quasimomentum, i.e., in the form of the commonly 
used dispersion law 



e( p) = — sin 7 tanh — sin - 

7 7 2 



2 P ■ , -2^0 

cos - + sinh — 

2 7 . 



1/2 



(12) 



A spinon corresponds in the usual Bethe ansatz classification of BAE solutions to a string of 
length 1 J7|. Naturally Eqs. (jl]) have string solutions of higher lengths too. Other spin ex- 
citations can be obtained as combinations of spinon quasiparticles and higher-length strings 
with different rapidities. However, spinons here are picked out because only their dressed 
energies may be negative, i.e., only spinons may form Dirac seas of the groundstate of the 
model. 

One can see that the dispersion law Eq. flljD of the low-lying excitation of the "easy- 
plane" two-chain spin ~ antiferromagnetic model is factorized into two parts: the gapless 
part at p = 0, tt and the gapful one at p — ir/2, cf. 0, |T0|| . The former corresponds to the 
oscillations of the magnetization, while the latter is connected with the oscillations of the 
staggered magnetization |J. The analysis, similar to the one of the solutions of Eqs. (0),(|D 
for nonzero magnetic field h 7^ (here we point out that according to the very accurate 
analysis |T(| the solution of the integral BAE in the first order approximation reproduces 
correctly both low- and high-coupling asymptotic behavior) shows that: (i) the dressed 
energy of a spinon as a function of the dressed quasimomentum has only one extremum, a 
maximum at p = tt/2 for 9 < 9 C and (ii) for 9 > 9 C there are two maxima and one minimum 
(situated at p = tt/2). At the (tri)critical point 9 C , the minimum disappears and two maxima 
joint into one more flat (at p — 7r/2). In the limit 9 —>■ 00 the mimimum is transformed 
into a cusp. It reveals that the gap of the staggered magnetization vanishes in this limit 
of two independent spin chains. This simple picture helps us to understand what happens 
if one switches on an external magnetic field h. Besides the usual phase transition to the 
ferromagnetic (spin-polarized) phase at 



(13) 



there is an additional transition between two phases. One of these corresponds to one Dirac 
sea of spinons (at small 9), while the other one is connected with two Dirac seas for the 
same kind of spinons (at large 9). It can also be seen from the r.h.s. of Eqs. (0), @ for 
the densities and dressed energies that the bare density and bare energy (corresponding to 
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terms which do not depend on p(u) and e(u)) have either one or two minima, respectively. 
Hence, they reproduce the same property in the dressed characteristics: The interaction 
simply "dresses" the (quasi)particles, as usual, but the "dressing" does not affect the picture 
qualitatively. The new critical field value can be approximated by h c « - sin 7 cosh -1 ^ in 
the first order approximation |J . In this approximation the tricritical point is the root of the 
equation 1 rs sinh ZL ^. At this point two second order phase transition lines h s and h c join. 
Hence, the "easy-plane" magnetic anisotropy in the antiferromagnetic two-chain model does 
not change qualitatively the groundstate behavior in the external magnetic field, cf. 0, |TU 



However it changes the critical values of the magnetic field and the intra-chain coupling. The 
difference between the two (gapless) phases is obvious: the first phase corresponds to the 
Neel-like antiferromagnetic groundstate for spins in both chains (along the zigzag line), while 
the second phase is connected with the Neel-like antiferromagnetic groundstates in each of 
chains, i.e. to effectively two magnetic sublattices in the two-chain model. 

That is why our simple model explains in which domains of parameters the two-chain 
spin system behaves like one-sublattice quantum "easy-plane" antiferromagnet, and where 
it behaves as the two-sublattice one. Note also that the phase transitions we study here are 
the manifestations of the commensurate-incommensurate phase transitions for spin systems. 
One can obviously see this, because the intra-chain coupling for two spin chains can be 
interpreted as the next-nearest neighbor spin interactions for a single spin chain of higher 
length N + + N_. Here the magnetic couplings are spin- frustrated, thus the emergence of the 
incommensurate magnetic states is understandable. 

As a consequence of the conformal invariance of (l+l)-dimensional quantum systems, 
the classification of universality classes is simple in terms of the central charge (conformal 
anomaly C) of the underlying Virasoro algebra fIT|| . The critical exponents in a conformally 



invariant theory are scaling dimensions of the operators within the quantum model. They 
can be calculated considering the finite-size (mesoscopic) corrections for the energies and 
quasimomenta of the groundstate and low-lying excited states. Conformal invariance for- 
mally requires all gapless excitations to have the same velocity (Lorentz invariance). The 
complete critical theory for systems with several gapless excitations with different Fermi ve- 
locities is usually given as a semidirect product of these independent Virasoro algebras. [ I 



Here we briefly sketch the procedure and write the results for the finite-size corrections to the 



energy, following the standard procedure, see, e.g. Refs. [flgl - One can see that for 9 < 9 C and 
for 9 > 9 C , h < h c , the conformal limit of our uniaxial two chain spin | model corresponds 
to one level-1 Kac-Moody algebra (one WZW model of level 1 with the conformal anomaly 
C = 1). The finite-size correction to the energy is rather standard (cf. 



E fs (N+ + AL) = ~v F + 2nv F (A l + A r ) , (14) 
o 

where v F is the Fermi velocity of the spinon and the conformal dimensions of primary 
operators are (please, pay attention: the lower indices denote the conformal dimensions for 
right- and left-moving quasiparticles, at the right and left Fermi point, respectively): 

/AM \ 2 
2A,, r = [ — ± zADj + 2n l>r , (15) 

where AM is an integer denoting the change of the number of particles induced by the 
primary operator, AD is an integer (half-integer) denoting the number of transfered particles 
from the right to the left Fermi point (backward scattering processes), n^ r are the numbers of 
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the particle-hole excitations of right- and left-movers. The values for the quantum numbers 
are restricted by AD = AM/2 (mod 1). The dressed charge z = £{Q) is the solution of the 
(standard) integral equation JO 



£(u) + / dvK(u - v)£(v) = 1 , (16) 
J (Q) 

taken at the limits of integration (these are the Fermi points, symmetric with respect to zero). 
In this phase there is only one region of integration over v. The dressed charge is a scalar. 
The behavior of our class of models in this phase in the conformal limit is rather standard 
[ PI . The correlation functions decay asymptotically oc (x — vpt)~ Al (x + v F t)~ Ar . The choise 
of the appropriate quantum numbers of excitations AM, AD and n^ r is determined for the 
leading asymptotics of correlators by taking the possible numbers with smallest exponents. 

But for 9 > 9 C , h > h c , the conformal limit of the "easy-plane" two-chain spin | model 
corresponds to the semidirect product of two level- 1 Kac-Moody algebras, both with confor- 
mal anomalies C = 1, i.e., to two WZW models both of level 1 || [TD[. The Dirac seas (i.e. 
the possible spinons with negative energies) are in the intervals [— Q + , — Q~] and [Q~,Q + ] 
(minima in the distributions of rapidities at =F#). This can be interpreted as the symmetri- 
cally distributed (around zero) Dirac seas of "particles" for [— Q + ,Q + ] and the Dirac sea of 
"holes" for [— Q~, Q~]. In fact the valley in the density distribution for "particles" and the 
maximum for "holes" are in the one-to-one correspondence with the maxima and minimum 
of the dispersion law for spinons. The second critical field h c in this language corresponds 
to the van Hove singularity of the empty band of "holes" . Naturally, the Fermy velocities of 
"particles" are positive, v F = (27ip(Q + ))~ 1 e'(u)\ u= Q+, while the Fermy velocities of "holes" 
are negative Vp = — (27Tp(Q~))~ 1 e r (u)\ u= Q-. The finite-size corrections to the energy for this 
case are 

E fs (N + + JV_) = -£(„+ + vp) + 2vr^+(A+ + A+) + v F (A T + A;)) , (17) 

where the dispersion laws of "particles" and "holes" are linearized about the Fermi points for 
each Dirac sea. The conformal dimensions of the primary operators are (the upper indices 
denote Dirac seas; the lower indices denote right and left Fermi points of each of these two 



Dirac seas, cf. |T0[ for the isotropic spin | two-chain model): 

(x_±AM+ - x + ±AM-) (z_±AD+ - z +± AD' 



2Af 

Lr 



2 det x 2 det z 



where the "minus" sign between the terms in square brackets corresponds to the right-, and 
"plus" sign to the left-movers. Here AM ± denote the differences between the numbers of 
particles excited in the Dirac seas of "particles" and "holes" , labeled by the upper indices. 
AD ± denote the numbers of backward scattering excitations, and the numbers of 

the particle-hole excitations for right- and left-movers of each of Dirac seas (for "particles" 
and "holes"). Please pay attention that AM^ and AD ± are not independent. Their values 
are restricted by the following connections: AM + - AM" = AM, and AD + — AD~ = AD, 
where AM and AD determine in a standard way the changes of the total magnetization 
and the total momentum of the system, respectively, due to excitations. Please note that in 
Refs. [IU|, these restrictions were missing; this resulted in, e.g., the invalid statement that 
four independent backscattering low-lying excitations are possible. However one can see that 
only two of them are really independent. The same is true for the excitations that change 
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the total magnetization of the system: there are only two independent of four possible such 
excitations. This is a direct consequence of the fact that only one magnetic field determines 
the filling of the Dirac seas for "particles" and "holes", or, in other words, two Dirac seas 
for spinons at ±9. 

The dressed charges Xik(Q k ) and Zik{Q k ) (i, k = +, — ) are matrices in this phase. They 
can be expressed by using the solution of the integral equation [^, [H| 

/(u|g ± )= (j Q \~j Q _)K(u-v)f(v\Q ± ) = K(u-Q ± ) , (19) 



with 18 



1 roo r—Q l 

z ik (Q k ) = S i>k + (-) k -( - )dvf(v\Q k ) 

2 JQi J~oo 

Xi k (Q k ) = 8 iik - R fc [ Q . dvf(v\Q k ) . (20) 



Notice, please, that the dressed charges depend on the value of the magnetic anisotropy 
7 via the kernels, while they depend indirectly on the value of the intra-chain coupling 
constant 9, only via the limits of integrations. In the first order approximation one can 
write the solutions as Xit{Q k ) ~ 8^ — (— ) k j® Q i dvK(v — Q k ) + . . . and z ik (Q k ) rs 8^ + 

(-)*(1/2)(Jq? -jZ^)dvK(u - Q k ) + .... The Dirac sea for "holes" disappears, naturally 
for h — > h c , 9 — > 9 C . The slopes of the dressed energies of "particles" and "holes" at Fermi 
points of the Dirac seas (Fermi velocities) differ in general from each other. Therefore we 
have a semidirect product of two algebras. Hence, in this region the dressed charges are 
2x2 matrices. This means that the conformal limit of the "easy-plane" two-chain spin 
| model corresponds to one or two WZW theories, depending on the values of the intra- 
chain coupling, magnetic anisotropy and magnetic field. At the critical line h c the Dirac 
sea of "holes" disappears as well as the components of the dressed charge matrix x (with 
square root singularities of the critical exponents for the correlation functions). Note that 
the dressed charge z becomes z = (2x)~ l at the phase transition line h c . This corresponds 
to the disappearence of one of the WZW CFTs. Unfortunately it is impossible to obtain an 
analytic solution to Eqs. ( |T9"D in closed form for a finite inter-chain coupling 9. Naturally 
in the limiting cases of two independent chains of lengths N±, 9 — > oo, and a single chain 
of length iV + + AL, 9 = 0, the solutions of Eqs. (|T^),(|l9|),(p0l) coincide with well-known 
ones, see Refs. ||18j| . The correlation functions of the uniaxial two-chain spin | model decay 



algebraically in this phase oc (x — Vpt)~ A i (x + v]?t)~ A i (x — Vpt)~ Ar (x + Vpt)~ Ar with the 
minimal exponents of possible quantum numbers of excitations AM ± , AD ± and nf r . We 
point out once more that the same magnetic field plays the role of a chemical potential for 
the "particles" and "holes" , or spinons of both Dirac seas in the second phase, and hence 
this choise of "minimal quantum numbers" is constrained. 

We must point out here that there is a crucial difference between our situation and the 
case of dressed charge matrices appearing for systems with the internal structure of bare 
particles [0. There the two Dirac seas of the groundstates are connected with different 
kinds of excitations, e.g., holons and spinons for the repulsive Hubbard model, or Cooper- 
like singlet pairs and spinons for the supersymmetric t — J model. They correspond to two 
different kinds of Lagrange multipliers, the chemical potentials and magnetic fields. Thus 
the low-lying excitations of the conformal theories in the spin and charge sectors of these 
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correlated electron models are practically independent of each other (spin-charge separation). 
Note that the spin and charge sectors are connected via the off-diagonal elements of the 
dressed charge matrix though. This is the consequence of the fact that say, holons or 
unbound electrons carry both charge and spin. On the other hand, two Dirac seas appear for 
the same kinds of particles for the models studied in this paper, which are also connected 
with the same magnetic field governing the filling of both Dirac seas. The latters appear 
due to two minima in the bare energy distribution and correspond to nonzero shift 9 in the 
Bethe ansatz equations. In other words, two Dirac seas are determined by the inter-chain 
coupling and appear if the values of coupling and external magnetic field are higher than 
the threshold values 8 C and h c , respectively. We believe that such a threshold behavior does 
not depend on the integrability of the model and is the generic feature for any multi-chain 
quantum spin models. 

The low temperature Sommerfeld approximation shows that as usually the low temper- 
ature specific heat is proportional to T out of critical lines. At the critical lines the van 
Hove singularities produce \JT low temperature behavior of the specific heat, while at the 
tricritical point we have T 1//4 behavior. 

What are the changes due to the different lengths of the chains AL 7^ AL? One can 
see obviously that the values of the momentum, energy and velocity of a spinon (which 
was v = (vr/7) sin 7 tanh(7r#/7)) become functions of iV + — AL. For example, the velocity 
renormalizes as v — > v[l + (AL — A^_) 2 tanh 2 (7r6 l /27)/A^ 2 ] _1 . This introduces dependences of 
the critical values 8 C and h c (as well as the saturation field h s ) on the difference AL — AL. 
Also, the Fermi velocities and Fermi points for the finite-size corrections become dependent 
on this difference. One can in principle consider different coupling constants J± for each of 
the chains (overall multipliers This produces the renormalizations similar to the action 

of AL 7^ AL, i.e., the velocity, e.g., renormalizes as v — » J+v[l + (J_/ J + ) 2 tanh 2 (7T0/27)] _1 . 



3 Connections to the Quantum Field Theories 

The studies presented in the previos section, being rather standard (note, though, some 
important new features, which were absent in the previous studies j|, |TT|, [IS], [| |T^, |T^] , such 
as the dependencies of the critical values of the inter-chain coupling and external magnetic 
field on the parameter of the magnetic anisotropy and on the difference in the lengths of 
the chains; also the important restrictions on the quantum numbers of low-lying conformal 
excitations). However we will use the results of that section for novel studies for a wider 
classes of exactly solvable models in Sections 3-5. For instance, in this section we point out 
the important similarities in the behaviors of the two-chain quantum spin model considered 
in the previous section and several models of QFT. 

Really, when examinating Eqs. fll]), one can see that these Bethe ansatz equations coincide 
with the ones, which describe the behavior of the spin (color) sector of some QFT. N± 
corresponds to the numbers of (bare) particles with the positive and negative chiralities. 
For example, for the chiral-invariant Gross-Neveu model |T3, 122]] we have to put 7^0 



(i.e. S{7(2)-symmetric case, equivalent to the S , f/(2)-symmetric Thirring model), and 6 = 
(1 — g 2 )/2g, where g is the coupling constant of the chiral invariant Gross-Neveu QFT Jl4 



As for the Lagrange multiplier h, it can play the roles of either an external magnetic field, 
or the chemical potential, or an external topological field, dual to the topological Noether 
current in QFT. Here we point out that in fact in QFT theorists are interested in physical 
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particles, which have the finite mass (gap). In the chiral-invariant Gross-Neveu model the 
gap of the staggered oscillations of the two-chain quantum spin model plays the role of 
the physical mass of the particle (spinor) [14], |I3| . As for the (gapless) oscillations of the 
magnetization of the two-chain spin model, we point out that they are the consequences of 
the lattice, and play the role of the massless fermion doublers of the lattice QFT p3 |. The 



results of the previous section mean that the behavior of the chiral-invariant Gross-Neveu 
model (or S'f/(2)-symmetric Thirring model) in an external magnetic field depends strongly 
on the coupling constant 9 (or equivalently on g). For 9 < 9 C the conformal limit of the 
QFT corresponds to one level- 1 WZW model with the conformal dimension C — 1. However 



for 9 > 9 C (-9 C - ^02 + 4 < 2g < -9 C + ^J9 2 c + 4) the conformal limit of this QFT in an 
external magnetic field corresponds to the semidirect product of two level-1 WZW model 
with the conformal dimensions C — 1. Two kinds of conformal points for this QFT have been 
mentioned already [p4] in a slightly different context. They were connected with one WZW 
theory or two WZW theories, coupled via a current-current interaction. This is related to 



right-left symmetry of the chiral invariant Gross-Neveu QFT (see, also, Refs. [32, |33] for the 
case of the QFT for the principal chiral field). 

Note, that the condition h > h c in the QFT means that the magnetic field is larger than 
the mass of the physical particle (color spinor). In this sense, in the region of magnetic field 
values h < h c the results of the QFT (see, e.g., [p2| ) predict zero magnetization, however 
a different lattice regularization, similar to the lattice scheme used in the previous section 
predicts a nonzero magnetization of the chiral-invariant Gross-Neveu model in this region. 
This is the indirect effect of the fermion doublers. In other words, it is connected with 
the well-known mapping of the lattice (e.g., Thirring) model under regularization on two 
continuum QFTs either both bosonic (free bosonic QFT and sine-Gordon one, f25|), or both 



fermionic ones (a free one and the continuum massive Thirring model). There are necessarily 
two such theories because of the Nielsen-Ninomiya fermion doublers: remember that we have 
started from a lattice [f23fl . 

For other models of QFT the procedure of the lattice regularization ^5J |27|, ^| has 
been used. Here 9 plays the role of the cut-off to preserve the mass of the physical particle 
to be finite. For example, for the [/(l)-symmetric Thirring QFT |^9[ ^] one can use the 
results of the previuos section with the limit 9 — > oo taken after the thermodynamic limit 
(L,N±,M — > oo with their ratios fixed, L is the size of the box). In this case one can 
obviously obtain the conformal limit of the theory with nonzero physical masses of the 
particles. Naturally in the limit 9 —>■ oo we ever exist in an external magnetic field in the 
phase with two Dirac seas. Here the latters correspond to the right- and left-moving particles 
(with positive and negative chiralities). Actually here our point of view coincides with the 
one of the field theorists. Recently it was shown in Ref. [30] that for (l+l)-dimensional sine- 
Gordon model the lattice regularization scheme in the "light-cone" approach gives similar 
to ours results for the conformal limit of the model. It was shown there that at the UV 
fixed point the conformal dimensions of the sine-Gordon model are determined by two U(l) 
charges of excitations (the usual one and the chiral charge). The chiral charge corresponds 
to the number of excitations transfered from one Dirac sea to the other, similar to our 
results (note that the above-mentioned lattice-regularized sine-Gordon case corresponds in 
our notations to 9 — > oo, where the integral equations for the particles with the positive 
and negative chiralities are totally decoupled). We point out here, that such a behavior is 
not unexpected, because the sine-Gordon QFT belongs to the same class of models, which 
is studied in our paper, i.e., its Bethe ansatz description features a shift of rapidities in the 
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Bethe ansatz equations in the lattice-regularized theory 



30 



4 Higher spin (chirality) generalizations 

For the higher spin generalizations of the Bethe ansatz theory presented in Section 2 we can 
write down BAE in the form 

M 

Y[e^(u Q ±6) = e™ M n up), (21) 

where n± = 2S± are the values of spins in each chain or the colors of bare particles in QFT. 
The eigenvalue of the transfer matrix can be writen as 



rr sinh(A - u a + + ^ M sinh(A ± 6) 

^ sinh(u a - A + 37-) ± sinl^ry^ - A =F 



X 



tt smh(u a - A + yy|) 
sinh(A - « Q - ryf ) 

Similar new phases with one or two kinds of Dirac seas for similar kinds of low-lying 
excitations exist also for a number of models in which n± ^ 1, e.g. for the higher-spin 
(S > |) two-chain models with equal spins in each chain, SU(n + 1) CIGN QFT |JT 



there n + — n_ = n ^ 1; for the principal chiral field models (nonlinear a-model) for CP- 
symmetric |32[ (there n + = n_ — > oo) and CP-asymmetric cases (there n + ^ n_, 
(n + + n_) — > oo, (n + — n_) fixed, i.e., the symmetry SU{2) x SU{2) oc 0(4)); and for the 
0(3)-symmetric nonlinear cr-model |34| as well as for spin-(S' + = 2n + ) - spin-(S'_ = 2n_) 
two-chain models {quantum two- chain ferrimagnet) . Note that for spins S ^ \ the procedure 
of the construction of the Hamiltonian is more complicated, because it corresponds to the 
two-chain uniaxial generalization of the Takhtajan-Babujian model, see e.g., Refs. ||35|| . For 



the simplest case of the isotropic exchange interaction between the spins and between the 
chains the Hamiltonian has the form: 



^ — X! ( ^ 2 (^5 + ,5+,ni,ni+l + 7^S_,S_ ,n 2 ,n 2 +l ,S—,ni,ri2 + 'Hs + ,S-,n 1 ,n 2 +l) + 

n ^ 

,£_,n 2 ,n 2 +l ,S-,ni,ri2 + 'Hs + ,S-,n 1 ,n 2 +l)} + 

,n 2 ,n 2 +l ,S— ,71i ,Tl2 + H S+ ,S_,ni,ra 2 +l )]) , (23) 

where [.,.] ({.,.}) denote (anti) commutator, 

■Si+S 2 j ^ S1+S2 x _ x 

n Sl ,S 2 ,n,n+l Yl S k 2 T m X IT ~ _ T > ( 24 ) 

3=|Si-S 2 |+l MS1-S2I+I «=|5i-5 2 | x i A/ 

x = Si t n02 t n+i, and 2xj = j(j + 1) — Si (Si + 1) — ^(^ + 1). The summation over n runs to 
N± in each chain. One can obviously see that for S± — = the Hamiltonian Eq. fl2"3] ) recovers 
the isotropic antiferromagnetic spin ~ Hamiltonian Eq. (|j) investigated in Section 2. For 
a single spin chain, 6 = 0, N + = N_ the Hamiltonian coincides with the known one of 



alternating spin chains [|36], |37], [38]. The Bethe ansatz studies of the model for n± can be 
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performed in the complete analogy with the above mentioned case n± = 1, keeping in mind, 
of course, the main difference: for the 5'[/(2)-symmetric or uniaxial higher spin models the 
groundstate corresponds to the filling up the Dirac seas for spin strings of lengths n±. ]35 



The well-known fusion scheme can be used for the case of a flavor-degenerate situation of 



the chiral invariant Gross- Neveu QFT, in the absence of flavor fields |39|]. Note that, except 
of the 0(3)-symmetric case, 7 = everywhere in the above-mentioned models of QFT. This 
corresponds to rational solutions of the Yang-Baxter equation for the two-particle scattering 
matrices. For the two spin chains the two-chain quantum ferrimagnet model corresponds 
to two Takhtajian-Babujian chains of different values of site spins coupled due to nonzero 
9. The total quasimomentum and the energy of the system in the framework of the lattice 



(local) regularization scheme for some QFT can be written as [23 



M Q 

-2a t E = ££ 7— N± In e n± (u a ± 6) 

± a =l aUa 
M 

iaP = EE iV ± ln e n± K± o) , (25) 

± a=l 

where a and a t denote the space and time lattice constants, respectively, and their ratio 
fixes the velocity of light ("light-cone" approach). The CP-symmetric (chiral invariant) 
case corresponds to the situation, in which n + = n_ = n. The Dirac seas are related to 
the dressed (quasi)particles with negative energies (strings of length n±). The behavior 
of the dispersion law for excited particle in the CP-symmetric case (n + = n_ = n and 
iV + = AL) is similar to Eq. fll2"l) : for instance, for the chiral-invariant Gross-Neveu QFT and 
principal chiral field model the r.h.s. of Eq. ( |1"2"D must be simply multiplied by sin(7rr/ri + 
1)/ sin(7r/n+l) and the parameter 6 in Eq. (|i~2j) has to be replaced by (n+ 1)0/2. r = 1, . . . , n 
is the rank of a fundamental representation of the SU(n + 1) algebra. All the previously 
mentioned characteristic features from the case n± = 1 persist. The differences are in 
the levels of Kac-Moody algebras in the conformal limit: The conformal anomalies are 
C = ^jg. Now the conformal field theory is a semidirect product of a Gaussian (C = 1) 
4(J and Z(n) parafermion models ||41|| : the operators identified from the scaling behavior of 



states consisting of Dirac sea strings only (found from finite-size corrections) are found to be 
composite operators formed by the product of a Gaussian-type operator and the operator 
in the parafermionic sector. To find a nonzero contributions from parafermions (constant 
shifts) one can consider the states with strings of other lengths then the Dirac sea present 
p2]j . For the scaling dimensions these shifts are ^+4 ? r = 1,2, . . .. 

From now on we concentrate on the n+ 7^ n_ situation. For the two-chain spin system 
the situation corresponds to the quantum ferrimagnet. Here we point out that due to the 
zigzag-like interactions in the system and spin frustration the ferrimagnets of this class are 
in the singlet groundstate (compensated phase) for h = 0, unlike the standard classical 
ferrimagnets in uncompensated phases. The integral equations that determine the physical 
vacuum of the systems are similar to Eqs. They reveal one or several minima of the 

corresponding distributions of dressed energies and densities with possible negative energy 
states, i.e., one or several Dirac seas: 

e T {u) + / dvK TT >{u - v)e T >{v) = h^-n T + ^ ^e°± 
J ± 

p T (u) + J dvK TT ,{u - v)p T ,{v) = ^Pr,± • (26) 
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The index r enumerates two possible Dirac seas and appears due to n + ^ n_, and the ± 
enumerate two possible minima due to the nonzero shift 9. The index r was naturally absent 
for the CP-symmetric case n + = n^. Note that for quantum two-chain ferrimagnets the 
investigated gapless phases in the groundstate in an external magnetic field are similar to 
the spin-compensated and uncompensated phases. Thus the phase transition between those 
phases is similar in nature to the well-known spin-flop phase transition in the classical theory 
of magnetism. Note, though, that the spin-flop transition is of the first order ("easy-axis" 
magnetic anisotropy), while the transition under study is the second-order one ("easy-plane" 
anisotropy). The Fourier transform of the kernel is given by 

2 coth{uj/2){diag{e- n+lu)/21 cosh(n + u;/2), e " n - |w/21 cosh(n_u;/2)) + 
a x {e- {n+ ~ n - )luj/21 - e -( n ++"-)l w / 2 l)] . (27) 

diag(a, b) is 2 x 2 diagonal matrix and a x is the usual Pauli matrix. Note, that after 
taking the limit (n + + n_) — > oo, which is the case of the CP-asymmetric case of the 
QFT for the principal chiral field, i.e. with the Wess-Zumino term the inverse kernel 
coincides formally (up to a constant multiplier) with the one for the case n + = n_ = 1. 
This indicates the global 0(4) (0(3)) symmetry of the principal chiral field |£J. There 



may be also two different behaviors, corresponding to one or several Dirac seas for n + ^ 1 
or n_ ^ 1. Naturally in the conformal limit the associated WZW CFTs have different 
conformal anomalies determined by n±: C± = ^2+2 ■ ^ or ^ ne determination of the Gaussian 
parts of the conformal dimensions of primary operators Eqs. ([18]) can be used. One has 
to add the input from the parafermionic sectors, too H, The elements of the dressed 
charge matrices are the solutions of the following system of integral equations: 



£t,t'{u) dvK T ,(ll - U)6r 1± (v) = &r,r' , (2£ 



in which the summation over ± is due to the two possible Dirac seas (two minima in the 
distribution of rapidities) at ±9. For different values of the spins, n + ^ n_, a transition 
between two different phases is induced by increasing an external magnetic field to some 



critical value, even in the absence of the shift 9 |3^, [3£|. It differs from the CP-symmetric 
case n + = n_, where the phase transition is only connected with the nonzero value of the 
intra-chain coupling parameter 9. For the CP-symmetric case one or two Dirac seas of the 
same type of excitations exist due to nonzero 9. But in the CP-asymmetric case the existence 
of two Dirac seas can be related to two kinds of different low- lying excitations (particles). 
They are strings of lengths n + and n_, respectively. In this situation the dispersion laws may 
be independent (not only factorized as for the previous CP-symmetric cases). The (new) 
phase transition at h c reveals the van Hove singularity of the empty Dirac sea for the longer 
strings. The spin saturation field h s is connected with the empty Dirac sea of strings of the 
smaller length. 



5 Multi-chain quantum spin models 

It is worthwile to mention that phase transitions in an external magnetic field, similar to 
the ones studied in this paper for uniaxial spin chains and QFT, have been already studied 



in the ID quantum alternating single spin chains |37, 38[ , spin 1 isotropic two-chain models 
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TOf , and correlated electron models with the finite concentration of magnetic impurities 
4^]. The Bethe ansatz equations for those models are similar to the ones studied in the 
present paper, Eqs. (TJ),(0). Note that the energies for spin models are defined (as usual for 
the lattice models) as first logarithmic derivatives of the transfer matrices. The factorization 
of the dispersion law for the lowest excitations (spinon) reveals essentially two kinds of 
magnetic oscillations: excitations of the magnetization and oscillations of the staggered 
magnetization, i.e., the manifestation of essentially two magnetic sublattices. Naturally, the 
existence of the latters persists in the continuum limit of such systems too (cf., for instance, 
with the standard theory of antiferromagnetism). Two non- ferromagnetic phases also reveal 
themselves in finite-size corrections to the energies of these quantum spin models. There 
instead of a scalar dressed charge for the phase with one Dirac sea for spinons, 2x2 dressed 
charge matrices appear in the second phase with two Dirac seas for for spin strings of different 
lengths in alternating spin chain |3"T| , 38 1, or for spinons of the same kind in zigzag-like coupled 
spin chains (see || [K| for the isotropic two-chain spin-1 model). 



The symmetry-breaking terms [the difference (n+— n_) = 2 (Si— S2), or nonzero 9] in BAE 
are actually the reason for the emergence of several gapless phases (or two Dirac seas) in the 
groundstate in an external magnetic field. It is also interesting to note that a homogenuous 
shift of rapidities can be removed for one Dirac sea for the periodic boundary conditions by 
an appropriate unitary (gauge) transformation (shift of variables), e.g., u a — > u a ± 9. But 
in the case of open boundary conditions, BAE take the form (for simlicity reasons we write 
the free boundary situation only, without the external boundary potential): 

Eel N ± (u a ±e) = Y[Y[e 2 (u a ±u,) . (29) 
± ± p 

It is clear that for the open chain one cannot remove the shift 9 of rapidities u a from one 
Dirac sea by a special choice of the gauge. From this point of view the latter case is close to 
the CP-asymmetric situation in QFT. 

One can see from the structure of the Hamiltonians that for the two-chain spin models 
the parameter 9 characterizes the intra-chain coupling for each chain (or the next-nearest- 
neighbor interaction in a single spin chain picture). It is obvious to introduce the series of 
{9j}J =1 (for each chain) and to construct the Hamiltonian of the exactly integrable multi- 
chain (J is the number of chains) spin model. For the simplest case of all S = | isotropic 
antiferromagnetic chains the Hamiltonian reads ||: 

PSn,pSn + l,p] + • • • + ( Ps n jS n J+l Ps n> jS n} J + l Ps n ,jSn+l,l J J J (^0) 

where A is the normalization constant (which depends on 9j), Ps a s b oc (1/2)1® I + 2S a <g> S& is 
the permutation operator and [., .] denotes a commutator. Note, that in the case of J ^ 2 the 
integrable model corresponds to the pair couplings not only between the nearest-neighbor 
spins but also to the next-nearest three spins, etc., couplings. All those terms are only 
essential in quantum mechanics, because in classical physics they are total time derivatives 



11] and do not change equations of motion. The Bethe ansatz equations have the form: 



J M 

JJ ef J (u m + 9 3 - 0i) = e mM J] e 2 (u m - u k ) , (31) 
j=i k 
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where M is the total number of down spins and Nj is the number of sites in the j-th chain. 
The previously studied situation J = 2 corresponds to the shift of the variables u m — > u m + 9 
with 9 2 — Q\ = —26. Now Oj — 1 determines the values of the intra-chain couplings in chain 

3- 

The analysis of the low-temperature thermodynamics of the multi-chain spin system is 
analogous to the situation of J = 2 studied in Sections 2-4. From the structure of the Bethe 
ansatz equations in the thermodynamic limit Nj,M — > oo, their ratios fixed, one can see that 
for J-chain model (for different Oj) there can exist, generally speaking, J phase transitions 
of the second order in the groundstate in an external magnetic field. These are nothing 
else than the commensurate-incommensurate phase transitions for the quantum multi-chain 
spin model with different couplings between the chains. The values of the critical fields 
h Cl , . . . , h Cj l and the value of the magnetic field of the transition to the ferromagnetic state 
h s depend on the set of Oj, i.e., on the intra-chain couplings (and also on the values of the 
magnetic anisotropy constants, which can be taken different for each chain; this does not 
destroy the integrability) . The ferromagnetic state is gapped, while all other phases are 
gapless in the integrable multi-chain spin quantum model. There are also J — 1 tricritical 
points at which the lines of the phase transitions h Cj join the line of the spin-saturation phase 
transition. Naturally, the phase that corresponds to the lowest value of the magnetic field, 
say h < h ci for special values of Oj (the condition is similar to < 9 C for J = 2), has in the 
conformal limit one scalar dressed charge. Hence, in the conformal limit our multi-chain spin 
model behaves as the level- 1 WZW CFT. In the next phase the multi-chain quantum spin 
model behaves as the semidirect product of two WZW CFTs, hence their dressed charges are 
2x2 matrices, and so on, until the last gapless phase, which corresponds to the semidirect 
product of J WZW CFTs with J x J dressed charge matrices. Note that J in this approach 
also denotes the number of possible Dirac seas (each of them is connected with the same 
magnetic field, so the excitations in each of them are not independent), and, thus, with 
one-half of the number of Fermi points. In the limit J — > oo (i.e. quasi-2D spin system) 
one obtains the (2D) Fermi surface instead of the set of ID Fermi points (the latters become 
disributed more closely to each other with the grouth of J). In this limit the differences 
between Oj tend to zero, and that is why the differences between Cj , h Cj and also between 
h Cj and h s disappear, too. Therefore in this limit the only h s survives. It means that for the 
quasi-2D limit of such an integrable model of J coupled quantum spin chains for J — > oo we 
expect only two phases in the groundstate in an external magnetic field: the ferromagnetic 
gapped one and the gapless phase, which in the conformal limit corresponds to one WZW 
CFT (with single scalar dressed charge). The phase transition between these two phases in 
the groundstate in an external magnetic field is of the second order. 

6 Behavior of the non-integrable multi-chain spin sys- 
tems 

So far we have studied only integrable multi-chain quantum spin models. We have shown 
that the commensurate-incommensurate phase transitions of the second order have to re- 
veal themselves in an external magnetic field due to the intra-chain interactions (or the 
next-nearest interactions in a single quantum spin chain picture). We have shown that the 
emergence of these phase transitions does not depend on the value of the site spins, they 
emerge in the presence of the "easy-plane" magnetic anisotropy, which keeps the system in 
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the critical (gapless) region. It is not clear, however, which features of the behavior of the in- 
tegrable models with the "fine-tuned" parameters have to exist for more realistic multi-chain 
models, and what are the qualitative differences, we expect to exist between the integrable 
multi-chain models and real multi-chain spin systems. 

We have to add one more thing to clarify the situation: We study (quasi) ID spin quantum 
models, for which one can use the Lieb-Schultz-Mattis theorem (and its generalizations) |^4, 
||. However, it is obvious that due to the frustration of the interactions between neighboring 
spins, and the presence of additional terms in the Hamiltonians, which violate the time- 
reversal and parity symmetries in the systems (spin chiralities or spin currents), for all spin 
models studied in the paper one cannot satisfy the conditions of the theorem. Hence it cannot 
be applied (at least directly). That is why for all the models we study there are no spin 
gaps (except for the trivial one for the spin-polarized groundstate). (Here we are not talking 
about the gaps connected with the magnetic anisotropy, but rather about the Haldane-like 
spin gaps |^5[ which appear even for the isotropic spin-spin interaction, and about fractional 
magnetization plateaux ||). As we argued before O], namely the presence of the chiral spin 
terms (or the operators of the nonzero spin currents) in the Hamiltonian (which are the total 
time derivatives and do not change the classical equations of motion but rather affect the 
topological properties, like the choise of the ^-vacuum in Haldane's approach) is the reason 
why the low-lying spin excitations (and particles for lattice QFT) for our class of models are 
gapless and our low energy theories are conformal. It has to be mentioned that recent results 
of the perturbative RG analysis of the zigzag spin ~ chain without three-spin terms shows the 
tendency the RG currents flow to the state with the parity and time-reversal violation 



By the way, one can obviously see that the XY limit of the two-chain spin model does not 
correspond to the free fermion point of the exactly solvable model, and this coincides with 
the results of Ref. [[SJ. Note, though, that in the latter it was erroneously concluded that the 



time-reversal and parity symmetries were violated by the two-chain zigzag spin Hamiltonian 
with only two-spin couplings (i.e. the nearest and next-nearest-neighbor interactions in the 
single chain picture), without spin current terms in the Hamiltonian. Hence the symmetry 
of the considered state was lower than the symmetry of the Hamiltonian there. 

Naturally, the relevant perturbations to our integrable models will immedeately produce 
spin gaps. As usual, the algebraic (power-law) decay of the correlation functions in the 
groundstate of the models considered in this paper determines the quantum criticality. This 
means that, starting from the (conformal) exact solutions obtained in this paper one can 
argue that the response of the more realistic spin systems to perturbations can be evaluated 
by using perturbative methods, e.g., in a renormalization group framework. For example, let 
us study the effect of relevant perturbations to the Hamiltonians considered, H r = H + 5 Hi, 
where one can choose as H r , e.g., the standard Heisenberg or uniaxial Hamiltonians for 
several coupled quantum spin chains, and as H the Hamiltonians of spin chains considered 
exactly in this paper for some values of 8, where the three-spin terms are relevant. The 
correction to the ground state energy and the excitation gap (mass of the particle in QFT) 
for the quantum critical system are: AE (x -S^ d+Z ^y, and m oc 8 l ' y i respectively, where d is 
the dimension of the system, and z is the dynamical critical exponent. For the conformally 
invariant systems studied here one has d = z = 1. The application of the standard scaling 
relations yields y+x = 2(= z+d), where x is the scaling dimension, i.e. x = 2A/ + 2A r , found 
in the previous sections (for the phases with the dressed charge matrices the summation over 
upper indices is meant). Hence the gap for the low-lying excitations (the mass of the physical 
particles in QFT) for the perturbed systems will be m oc S l / 2 ( l ~ A i-^ r ) . Note that because 
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of scaling, the behavior of the critical exponents (which are related to the exponents we 
introduced for the integrable multi-chain spin models) in the vicinities of the lines of the 
phase transitions has to be universal, and this can be checked experimentally. We expect 
that the spin gap has to exist for the values of the isotropic zigzag inter-chain coupling 
higher or of order of 0.5 for the two-chain spin | system |J, where the three-spin couplings 
are relevant and the emergence of the spin gap is known exactly [47]. 

Very recently, the density matrix renormalization group numerical studies of the two- 
chain zigzag spin | model (without chiral three-spin terms in the Hamiltonian) were per- 
formed [Q. These numerical studies strongly support the picture proposed here (see also 
Ref. f9j): the magnetization as function of the magnetic field in the groundstate reveals (i) 
one second order phase transition (to the spin-saturation phase) for the weak intra-chain 
coupling; (ii) one more second order phase transition between the magnetic (gapless) phases 
in the intermediate region of the intra-chain coupling and (iii) in addition to those second 
order phase transitions, one to the gapful phase with zero magnetization (plateau) for the 
intra-chain coupling value of 0.5. 

We should also mention that it is not the chiral spin terms (as implied in Ref. [jl0[j ) but 



the intra-chain coupling that is responsible for the commensurate-incommensurate phase 
transitions between the gapless phases in this class of models. As for the aforementioned 
spin currents, their "fine-tuned" values produce the cancellation of the spin gap for zero 
magnetic field ||. We should also note that to our mind some features of the phase diagram 
obtained in Ref. [pl| are artifacts of the small number of sites involved into the numerical 
calculations. In Fig. 5 of Ref. |TP|] in the regions of 0.52 < k, < 0.6 (corresponding to 
intra-chain couplings, normalized to the value of the inter-chain interaction, in the domain 
[0.54-0.75]) we can obviously see that when increasing the value of the magnetic field one 
goes from the gapped phase with zero magnetization into the gapless one with two Dirac seas 
of the low-lying excitations, then reaches the gapless phase with one Dirac sea, then returnes 
to the gapless phase with two Dirac seas, and finally reaches the spin-saturated phase. To 
our mind this return to the already passed phase is non-physical. One can clearly see that 
the region for the values of the intra-chain couplings, where these strange returns happen, is 
reduced when going from 16 sites in numerical calculations to 20 sites. This confirms that 
presently achieved sizes of the quantum systems for numerical calculations can produce even 
qualitatively invalid results, and analytic calculations are necessary, too. 

We point out that despite the fact that the relevant perturbations in general produce a 
gap for the low-lying excitations, one can apply the results of this paper to the real gapless 
multy-chain spin systems, too. For example, it was recently observed that even for the 



two-leg ladder system SrCai 2 Cu 24 04i the spin gap collapses under pressure. |49 



7 Concluding remarks 

In this paper, motivated by recent progress in the experimental measurements for multi- 
chain spin systems, we have theoretically studied the behavior in an external field of a wide 
class of the multi-chain quantum spin models and quantum field theories. First, we have 
investigated the external field behavior of the exactly integrable two-chain spin | model 
and have shown that the inclucion of the magnetic anisotropy of the "easy-plane" type, 
with which the system stays in the quantum critical region, does not qualitatively change 
the behavior in an external magnetic field. However, we have shown that the magnetic 
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anisotropy changes the critical values of the magnetic fields and intra-chain couplings, at 
which the phase transitions occur, and affects the critical exponents. We have shown that the 
external-field-induced phase transitions we discussed are the commensurate-incommensurate 
phase transitions due to the next-nearest-neighbor two-spin interactions, which are present in 
these multi-chain models with zigzag-like couplings. We have pointed out that the low-lying 
excitations of the conformal limit of our class of multi-chain spin models are not independent 
in the incommensurate phase, because they are governed by the same magnetic field. We 
have shown that these two-chain quantum spin models share the most important features of 
the behavior in an external field with the wide class of the (1+1) quantum field theories. 

We have introduced higher-spin versions of the two-chain exactly solvable spin models, 
e.g., we have investigated the important class of ID two-chain quantum ferrimagnets with 
different spin values in the sites of each chain. Here we have shown that the phase transitions 
in an external magnetic field in this exactly solvable two-chain quantum ferrimagnet are 
similar in nature to the phase transitions between the spin-compensated and uncompensated 
phases in ordinary classical 3D ferrimagnets. 

We have also studied the behavior of the multi-chain exactly solvable spin models in 
an external magnetic field, and shown how the additional phase transitions arising due 
to the increasing number of chains vanish in the quasi-2D limit. Hence, to the best of 
our knowledge, we have proposed the first exact scenario of the transition from ID to 2D 
quantum spin models in the presence of an external magnetic field. We have argued that the 
commensurate-incommensurate phase transitions in the multi-chain quantum spin models 
have to disappear in the limit of an infinite number of chains. 

Finally, we have shown how the relevant deviations from the integrability, e.g., the ab- 
sence of the three-spin (spin chiral) terms in the Hamiltonians, which separately break the 
parity and time-reversal symmetries, give rise to gaps in spectra of the low-lying excita- 
tions of the multi-chain quantum spin systems and we have calculated the critical scaling 
exponents for these gaps. We pointed out the qualitative agreement of our exact analytic 
calculations with recent numerical simulations for zigzag spin models. 

I am grateful to A. G. Izergin, S. V. Ketov , A. Kliimper, V. E. Korepin, G. I. Japaridze, 
A. Luther and A. A. Nersesyan for helpful discussions. I thank J. Gruneberg for his kind 
help. The financial support of the Deutsche Forschungsgemeinschaft and Swedish Institute 
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